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X-ray reflection in a sample of X-ray bright Ultraluminous X-ray 
sources 
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ABSTRACT 

We apply a reflection-based model to the best available XMM-Newton spectra of X-ray bright 
UltraLuminous X-ray (ULX) sources (NGC 1313 X-1, NGC 1313 X-2, M 81 X-6, Holm- 
berg IX X-1, NGC 5408 X-1 and Holmberg II X-1). A spectral drop is apparent in the data 
of all the sources at energies 6-7 keV. The drop is interpreted here in terms of relativistically- 
blurred ionized reflection from the accretion disk. A soft-excess is also detected from these 
sources (as usually found in the spectra of AGN), with emission from O K and Fe L, in the 
case of NGC 5408 X-1 and Holmberg II X-1, which can be understood as features arising 
from reflection of the disk. Remarkably, ionized disk reflection and the associated powerlaw 
continuum provide a good description of the broad-band spectrum, including the soft-excess. 
There is no requirement for thermal emission from the inner disk in the description of the 
spectra. The black holes of these systems must then be highly spinning, with a spin close to 
the maximum rate of a maximal spinning black hole. The results require the action of strong 
light bending in these sources. We suggest that they could be strongly accreting black holes in 
which most of the energy is extracted from the flow magnetically and released above the disc 
thereby avoiding the conventional Eddington limit. 

Key words: Une: formation - black hole physics - X-rays: galaxies - X-rays: general 



1 INTRODUCTION 

Ultra-luminous X-ray (ULX) sources are point like, nonnuclear 
sources observed at other galaxies, with observed luminosities 
greater than the Eddington luminosity for a lOMp st ellar mass 
black hole (BH), with Lx^lO^' ergs"' ('Fabbiano"l989^. The true 
natur e of these objects is still open to debate ( Miller & Colbert! 
|2004) . One fundamental issue is whether the emission is isotropic 
or beamed along our line-of-sight. A possible scenario for geomet- 
rical beaming involves super-Eddin gton accre tion during phases of 
thermal-timescale mass transfer (Kinal l2002h . Alternatively, if the 
emission is isotropic and the Eddington limit is not violated, ULX 
must be fuelled by accretion onto Intermedi ate-Mass BH (IMBH), 
with masses in the range 100-10 000 Mq jColbert & Mushotzkvl 
Il999h . Currently, there is no agreement regarding the nature of 
these sources. It is possible that some ULX appear very luminous 
because of a combination of moderately high mass, mild beam- 
ing and mild super-Eddington emission. It is also possible that 
ULX are an inhomogeneous population, comprised of both a sub- 
sample of IMBH and moderately beamed stellar mass black holes 
jpabbiano & Whitdl2 006: Miller & Colb erlll2004h . 

NGC 5408 X-1, HOLM II X-1, M 81 X-9 (which is in the 
companion-dwarf galaxy Holmberg IX, hereafter called HOLM IX 
X-1), the ULX in NGC 1313 and M 81 X-6 are ULX located 
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in dwarf (NGC 5408 X-1, HOLM II X-1 and HOLM IX X- 
1) and spiral (NGC 1313 X-1, NGC 1313 X-2 and M 81 X-6) 
galaxies, respectively. All these ULX peak in X-ray luminosity 
above Lx = IxlO''" ergs"', thus being excellent targets for test- 
ing spectral models to the data. They are n earby and located at 
distances of D = 4.8, 3.5, 3.7, 3.63, 4.50 Mpc jKarachentsevll2002l : 
IStobbart et al]|2006l : lLiu & Pi Stefandl2008l) . for NGC 5408 X-1, 
HOLM IX X-1, NGC 1313 X-1, M 81 X-6 and HOLM II X- 
1 respectively. NGC 5408 X-1 is among the few ULX for which 
(10-200 mHz) Quasi-Per i odic Oscillations (QPO ) have been found 
dStrohmaver et alj|200l IPewangan et al.ll2006h . NGC 5408 X-1 
exhibits X-ray timing and spectral properties analogous to those ex- 
hibited by Galacti c stellar-mass black hole i n the very high or steep 
power-law state ( iRenullard & McClintockl 1 200 6 i. but with the 
characteristic variability t imescales (QPO and br eak frequencies) 
consistently scaled down jStrohmaver et al.ll2007h . For NGC 5408 
X-1 the inferred characteristic size for the X-ray emitting region is 
» 100 times larger than the typical inner disk radii of Galactic black 
holes. The highest signal-to-noise spectra of ULX can often be fit- 
ted by composite models of a thermal disk together with a hard- 
powerlaw tail, with a low measured disk temperature of »;0.2keV 
(iMilleret alJ2003ll2004h . If such factors are entirely due to a higher 
black hole mass, they imply an IMBH with MalOO - 10000 Mq. 
Both the morphology and flux of the optical high-excitation neb- 
ulae detected around some ULX probabl y rule out strong beam- 



ing as the origin of the X-ray emission jPakull & Mirioni 20o3; 
iKaaret et al.ll2004l : ISoria et alfcood i FAbolmasov et al.ll2007n . 
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The detection for most ULX of spectral c urvature, in the form 
of a deficit of photons at energies E^lks W ( [Roberts et al.ll20o3 : 
IStobbart et"aDl2006l : iMiv awaki et al.' 2009) has led to the sugges- 
tion that most ULX have spectral properties that do not correspond 
to any of the accretion states known in Black Hole Binaries (BHB), 
making it unlike ly that ULX are pow ered by sub-Eddington flows 
onto an IMBH ^Roberts et alj 2007[). The application of Comp- 
toniz ation models to the data (Stob bart et alj2003 : lGladstone et all 
1200^ and references therein) results in strikingly high and low val- 
ues for the coronal opacity (t>5) and the electron temperature 
{kTc =1-3 keV), difficult to explain for the expected physical 
conditions in a corona surrounding the black hole. This is very dif- 
ferent to the typical values found for BHB during the low/hard state, 
with spectra dominated by ComptonizationQ and appears irrecon- 
cilable with the IMBH model, which assumes that they operate as 
simple scaled-up BHB. 

Here, we present an alternative interpretation of the spec- 
tral shape, based on a physically-justified model commonly used 
on other accreting black holes. The soft part of the spectrum 
(£<1 keV) - the soft X-ray excess - and the high-energy curvature 
are j ust aspects of a reflection spectrum expected from accretion 
( IGui lbert & Rees 1988; George & Fabian 1991). A major compo- 
nent of reflectio n, the broad iron K line, has been foun d in many 
Seyfert galaxies dTanaka et al.lll995l.lNandra et al.ll2008h. a ccreting 



stellar-mass black holes jMiller et al 20071: Reis et alj|200 9a) and 



even accreting neutron stars teackett et al. l2008l : lReis et ali2009bi) 



Both the soft excess and the relativistic broad iron K line have 
recently been demonstrated to be part of the same physical pro- 
cess, i.e., the reaction of the disk to irradiation from a high-energ y 
source, in the Seyfert-1 galaxy IH 0707-495 dFabian et alj|200^ . 
In this paper we investigate whether reflection models can account 
for the spectra of these ULX and, if so, we put them in the context 
of the accreting black holes known so far. 

For this study, we choose the ULX with the best available data 
and with the longest exposure time observations (^lOOks of expo- 
sure time) of the XMM-Newton satellite. By using only the highest 
quality data from the widest band pass, highest sensitivity instru- 
ments available we expect to make good statements on the accretion 
processes in these ULX. In Section[2]we describe the observations 
and data used, in Section[3]we report on the results of the applica- 
tion of fits with the reflection model to the spectra and in Section|4] 
discuss the results obtained. 



2 OBSERVATIONS AND DATA REDUCTION 

In this work, we consider the time averaged XMM-Newton EPIC- 
pn spectra from the longest available observations of sources with 
a rate >0.5 counts/s in the X-ray band. The datasets were obtained 
through the XMM-Newton public data archive. The EPIC-pn cam- 
era has a higher effective area than the EPIC-MOS cameras, and 
drives the results of any joint spectral analysis. The reduction and 
analysis reported in this work used SAS version 8.0.0. We checked 
for pile-up in all the observations and found that this was not sig- 
nificant (i.e. less than 5 per cent for the high energy channels) for 
all the observations. 

In Table [T] we present a log of the observations. We applied 

' It has to be noted that similar values, i.e. kTg = 1 - 3keV and r>5, 
have been found in the description of the spectra of BHB durin g the steep 
powerlaw state (e.g. GRO J1655-40 and GRS 1915+105: Makishima et alj 
l2000tlKubota & Donel200llUeda et al.l2009t) . 



the standard time and flare filtering (rejecting high-background pe- 
riods of rate >0.4 counts/s, as recommended for the pn camera|3). 
We filtered the event files, selecting only the best-calibrated events 
(pattern<4 for the pn), and rejecting flagged events (flag= 0). 

For each exposure, we extracted the flux from a circular re- 
gion on each source. The background was extracted from a circular 
region, not far from the sources. Both the coordinates of the cen- 
troids and th e radius used for the re gions were the same as those 
employed bv lGladstone et al.l ( |2009|) . We built response functions 
with the SAS tasks rmfgen and arfgen. We fitted the background- 
subtracted spectra with standard models in XSPEC 12.5.0 lAmaudl 
[1996). All errors quoted in this work are 90% confidence errors, ob- 
tained by allowing all variable parameters to float during the error 
scan. Owing to the uncertainties in the EPIC calibration, we used 
only the 0.3-10 keV range. The resulting spectra were grouped with 
the FTOOL grppha to bins with a minimum of 20 counts each. 



3 SPECTRAL ANALYSIS AND RESULTS 

In contrast to previous work, where Comptonization models are 
used to fit the data, we looked at the observations of NGC 1313 
X-1, NGC 1313 X-2, M81 X-6, HOLM IX X-I, NGC 5408 X-I 
and HOLM II X-1 from the point of view of a broad iron emission 
line, i.e. a reflection spectrum. An absorbed powerlaw was fitted to 
the data in the bands 1.2-3 and 8~10keV. In the case of NGC 1313 
X-2 and M81 X-6 we used the energy range 1-2 and 8-10 keV 
instead. The ratio of the whole 0.3-10 keV data set to the absorbed 
powerlaw is shown in Figure [T] The column densities found are in 
the range of A'h = (0. 1 - 0.5)x 10"^ cm"^, thus lar ger than the corre- 
sponding galactic Hydrogen column densities jpickev & Lockmanl 
I1990i) . indicating the presence of either local extra-absorption to 
the sources or intrinsic to the host galaxies. A large s kewed emis- 
sion feature, similar to those found in BHB and AGN jMiller et al.l 
^007; Tanaka et al. 1995; Nandra et al. 2008) is evident about 3 and 
8 keV, and a soft excess (like those found in AGN) below 1 keV. 

The data between 1.5 and lOkeV wer e then mode lled with 
a powerlaw and a relativistic emission line ( lLaoilll99lh . A good 
fit is obtained with F = 1.5-3 (see Table O. The addition of 
this line represents an improvement in the statistics of A;^^ = 
60, 115,367, 103,65, 10 (for 5 extra d.o.f.) for NGC 1313 X-1, 
NGC 1313 X-2, M81 X-6, HOLM IX X-I, NGC 5408 X-1 and 
HOLM II X-I, respectively. The line was thereafter significant in 
all the sources with exception of HOLM II X-1. HOLM II X-1 has 
the lowest laor line detection, due to both the short exposure time 
and steepness of the spectrum (we will refer to this point later). 
Further longer exposure observations of HOLM II X-1 are needed 
to fully address this point. Notice that no extreme inclination is re- 
quired to reproduce the spectral drop (r;40 - 60°). The strength of 
the line (equivalent width), the rest frame energy and the emissiv- 
ity index of the line are in the range of EW = 0.15 - 1.27 keV, 
£ = 6.4 - 6.97 keV and 15- = 6 - 10 (see Table[2). For the case of 
NGC 1313 X-1, NGC 1313 X-2, M81 X-6 and HOLM IX X-1 the 
rest-frame energies imply that we are dealing with ionized reflec- 
tors (the steepness of the powerlaw for the case of NGC 5408 X-1 
and HOLM II X-1 did not allow us to obtain reliable values for 
the line parameters). There are remarkable emission line spectral 
features around £a0.6 and 1 keV in the spectra of NGC 5408 X-1 
and HOLM II X-1 (see Figure[T), which seem likely to come from 
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Table 1. XMM-Newton observations log and effective exposure times of XMM-Newton instruments. 



Source 


ObsID 


Date 


Exposure time ( s) 


Rate (counts s ' ) 


NGC 1313 X-1 


0405090101 


2006-10-16 


121212 


0.73 


NGC 1313 X-2 


0405090101 


2006-10-16 


121212 


0.67 


HOLM IX X-1 


0200980101 


2004-09-26 


117 023 


1.54 


NGC 5408 X-1 


0302900101 


2006-01-15 


130549 


0.96 


M81 X-6 


0111800101 


2001-04-23 


132 663 


0.48 


HOLM II X-1 


0200470101 


2004-04-15 


104 677 


3.13 



VIII K and Fe XVII and XVIII L. These lines have been usu- 
ally understood as coming from the hot diffuse gas around these 
sources and/or emission from a photoionized stellar wind, similar 
to what is seen in some high-mass X-r ay binar ies in our own Gal axy 
jpeng & KaaretlbOOAlstobbart et alJ[2006l) . iGoad et al.l j2006l) ar- 
gued that these features commonly seen in the spectra from ULX 
correspond to poorly modeled spectra of ULX considering ab- 
sorbers with solar metal abundances (i.e. they instead found that the 
RGS spectrum of HOLM II X-1 is consistent with an absorber of 
0.6 times solar abundance). Nevertheless, this value does not mean 
anything about the local environment of the source, since this is 
the metal abundance found in the total gal actic column density i n 
the direction to the source As pointed bv'R oss & FabianI ( |200^ . 
these lines can plausibly come from the disk of the sources. They 
originate through X-ray ionized reflection, which occurs when a 
surface is irradiated with X-rays. The reflection continuum is then 
complex and needs to be modelled properly. 

3.1 Reflection spectra 

We have therefore fitted IRoss & Fab iai? ionized reflection 

models to the whole 0.3-10 keV observable band. Very good agree- 
ment is obtained with a model which is relativistically blurred, 
has supersolar abundances between 4 and 10 and photon index of 
r«1.5 - 1.6 for NGC 1313 X-1, NGC 1313 X-2 and HOLM IX 
X-1, r«2 for M81 X-6 and steeper (r>2.5) for NGC 5408 X- 

1 and HOLM II X-1 (see Table |2). In the case of NGC 5408 
X-1 and HOLM II X-1, this model has worse statistics ix' Iv = 
736/577,945/779), presumably due to the steepness of the spec- 
trum, which makes the statistics be dominated by the softest part of 
the spectrum (£'<1.5keV). This would make easy the determina- 
tion of the O abundance (because the lowest part of the spectrum is 
very well sampled) at the expense of having Fe badly constrained 
(being responsible for the emission at the highest part of the spec- 
trum). Nevertheless, a study with variable abundances is beyond the 
scope of this paper. This model provided a good description of the 
data (see Table|2l. The ionization parameter of the reflection model 
is in the range of f = 29 - 2 500. The high values found for the 
ionization parameter for NGC 1313 X-I, NGC 1313 X-2, M 81 
X-6 and HOLM IX X-I are very likely related to the featureless 
low energy (£< 1 keV) spectra for these sources. 

The relativ istic blurring is accomplished using a Kerr ker- 
nel ( lLaoJll99lh with a radial emissivity index in the range of 
q = 4-9 (where the emissivity has the radial dependence r^'') from 
Rin = 1.2- lARg for NGC 1313 X-1, NGC 1313 X~2 andM81 X- 
6 and = 2.0 - 2.5 Rg for HOLM IX X-1, NGC 5408 X-1 and 

^ In this paper we apply a model describing the X-ray emission local to 
the source which incorporates the metal abundances as a parameter, thus 
potentially allow us to unveil the metallicities local to these sources. 



HOLM II X-1. Here, iRg = GM/c^. The outermost radius is fixed 
to 400 a canonical value inferred from fits to spectra of AGNs. 
Since there is no previous measurement for the inner disk inclina- 
tion so far, it was constrained in the fits to be within the range of 
/ = 40 - 60°, i.e. an intermediate value. The inner disk radii found 
for the former ones indicate that the matter reaches the last stable 
orbit for a very rapidly spinning black hole, as inferred from the 
highly skewed profiles for the broad Fe K line for these sources (see 
Figure [T]l. They show very strong reflection components as well 
(^icfl/^totai = 0.7 - 1), i.e. the spectrum is dominated by reflection, 
thus being reflection-dominated sources. Indeed, NGC 1313 X-2 
and M81 X-6 are the most extreme cases, with no signs of emis- 
sion from a powerlaw high-energy source. In contrast, for the cases 
of HOLM IX X-1, NGC 5408 X-1 and HOLM II X-1, the spec- 
tra are powerlaw dominated (powerlaw -dominated sources). Curi- 
ously, the reflection dominated sources are the ones which show 
the smallest inner disk radii, perhaps indicating that light bending 
is an important effect for them. It is worth noting that for the most 
reflection-dominated sources, i.e. NGC 1313 X-2 and M81 X-6, 
the residuals at the broad Fe complex look peaky and that this is 
a natural consequence of the fact that we are modelling with a 
uniform disk, which could differ in properties to its outer region 

0. In contrast, for the powerlaw dominated sources, the inner disk 
radii are larger In the case of HOLM IX X-1, the larger radius 
{Rm = 2.0 Rg) is accompanied by the flattest powerlaw in the over- 
all sample (F = 1.504±0.004), and this could imply that reflection 
occurs farther from the black hole with respect to the reflection 
dominated ones. In the case of NGC 5408 X-1 and HOLM II X- 

1, with steeper powerlaw indices (r>2.5) and showing worse fit 
statistics, due to the lower signal-to-noise of the high-energy chan- 
nels, no definite statements can be made. In none of the spectra 
does the addition of a disk component (diskbb in XSPEC) im- 
prove the statistics. Therefore, we conclude that we cannot detect 
thermal emission from the disk for these sources (we shall return to 
this point in Section|4j. 

The fits with both the laor and the reflection models above 
indicate high values of the emissivity parameter (q = 4 - 9), indi- 
cating that relativistic blurring is an important effect and the emis- 
sion is concentrated near the innermost radii. The respective in- 
ner disk radii obtained are low (/?,„ = 1.25 - 2.5 Rg), indicating 
X-ray emission from regions very close to the black hole, where 
strong light bending dominates. These effects mean that high val- 
ues of the spins of the black holes are expected. To study this pos- 
sibility further, we substituted the Kerr kernel (kdblur), which 
assumes maximal spin, by the relativistic kernel of (kerrconv; 



The emissivity parameter is expected to vary between the inner and the 
outer boundaries of the disk, giving very smooth (high y) and peaky (low 
q) components, respectively. In this paper we just model the smooth com- 
ponent, expected to come from the inner region of the disk. 
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NGC 1313 X-1 



NGC 1313 X-2 




Energy (keV) 



Energy (keV) 



Figure 1. Ratio of the XMM-Newton pn spectra of NGC 1 3 1 3 X-1 , NGC 1313 X-2, M8 1 X-6, HOLM IX X-1 , NGC 5408 X-1 and HOLM 11 X-1 (upper-left 
to lower-right) to the powerlaw which fits best between 1.2-3 and 8-10 keV. A broad iron line and a soft excess are evident. Data have been rebinned for visual 
clarity. 



iBrenneman & Reynold j 1200^ , which assumes variable spin pa- 
rameter (a - cJ/GM^; where J is the angular momentum of a black 
hole of mass M). We assumed the inner disk radii to be equal to the 
marginal stability radius (i.e. the closest inner disk radius to the 
black hole) and constrained the disk inclination to an intermediate 
value ((■ = 40 - 60°). Thanks to the flat spectra of NGC 1313 X-1, 
NGC 1313 X-2, M 81 X-6 and HOLM IX X-1, which lead to a 
strong signal at high energies, we could constrain well the values of 
the emissivity parameter of the disk and the spin of the black hole. 
The spins for the black holes of these sources are very high, and 



close to the value of a maximal spinning black-hole (a = 0.998), 
with M 81 X-6 hosting the most rapidly spinning black hole of 
the sample (a = 0.993 - 0.998). In the case of HOLM II X-1, the 
lowest value found (a = 0.47+0.12) might be a consequence of the 
steep powerlaw and short exposure, which lead to little high-energy 
signal (NGC 5408 X-1 is similar). Future work is needed in order 
to fully address this point. As in the application of the Kerr ker- 
nel, the values obtained for the ionization parameters are higher for 
NGC 1313 X-1, NGC 1313 X-2, M 81 X-6, HOLM IX X-1 than 
in NGC 5408 X-1, HOLM II X-1, which correspond to featureless 
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low-energy spectra for the former ones due to the current exposure 
times used. The opposite occurs for NGC 5408 X-1 and HOLM II 
X-1, for which we could constrain very well the ionization param- 
eter (i.e. the disk density) to ^ = 49.4±0.8, 29±6, due to the good 
signal of the spectrum at low energies (£<1.5keV). This value is 
low for both sources and indicates a high density for the reflecting 
medium (i.e. the disk). 

For all the sources, we find that iron is overabundant com- 
pared to the Sun by 4 - 9 times. Usually, ULX have been asso- 
ciated with metal poor environments, because stars in these sites 
lose less mass during their evolution and can reach very high 
masses and collaps e to high-mass black h oles (M<80Mo) at the 
end of their lives ( Fryer & Kalogerj|2001 1. Furthermore, there is 
observational evidence of real low-metal abundances in the galax- 
ies which host these ULX. Metal abundance measurements from 
the HII regions in NGC 1313 and HO LM II X- 1 have values 
Z^ Q QOS, OA fa adfield & Crowther 2007 dPakull & Mirioni 2002) 
and lMapelli et al.l ^2009) have found similar values for other ULX 
host galaxies. These facts agree with the known mass-metallicity 
relation, in which t he smallest galaxies are the poorest in metals. 
IWinter et al.l j2007h found metal overabundance for a sample of 
ULX, with mean metallicity close to the solar value 0. With ex- 
ception of HOLM II X-1 and HOLM IX X-I, where the metal 
abundance they obtain is compatible with ours, we can explain the 
discrepancy in the level of metals by means of the different model 
used for the soft excess (thermal emission from an accretion disk 
in their case). Nevertheless, they found, like in our study, that the 
galaxies do not obey the mass-metallicity relation in X-rays. This 
does not mean a conflict with optical measurements, since we are 
measuring metallicities local to the sources in X-rays. One possi- 
ble scenario to explain such high values for the metallicities ob- 
served local to the sources, without any apparent chemical metal 
enrichment in their surrounding environments, is that the metals 
were blown out by the type II SNe event and then recaptured (sim- 
ilar to the fallback disk model of .Li.2003) . We shall return to this 
point later in Section]?] 



4 DISCUSSION 

Remarkably, a relativistically-blurred, reflection-dominated model 
gives a very good fit to the sample of ULX with best quality ( 1 00 ks 
of time exposure) XMM-Newton data. Very skewed iron line pro- 
files have been found, implying that the emission region is very 
close to the black-hole and this fact allowed us to determine their 
spin. We find that all the ULX in our sample are close to maxi- 
mally rotating, with the exception of NGC 5408 X-I and HOLM II 
X-1, for which the steep powerlaw did not allow us to properly 
determine parameters from the Fe line. We have found that, the 
sources showing a reflection dominated spectrum (NGC 1313 X-I, 
NGC 1313 X-2 and M 81 X-6) have maximally spinning black- 
holes. This fact can be understood if light bending is an important 
effect for these sources. For these sources, the high-energy source 
(i.e. the jet or corona) is very close to the black hole and the ob- 
served direct continuum (i. e. powerlaw) is very low . This would 
correspond to Regime I of iMiniutti & FabiMil ( |2004[) . correspond- 
ing to a low height of the primary source, where strong light suf- 
fered by the primary radiation dramatically reduces the observed 



^ The sources with the longest exposure times in their observations, i.e. 
HOLM II X-1 and HOLM IX X-1 clearly are outliers of this relation. 



powerlaw emission component at infinity. HOLM IX X-I shows 
a powerlaw dominated spectrum with a slightl y narrower Fe line 
profile. This would correspond to Regime II of'Miniutti & Fabiai] 
((2004), when the height of the primary source is larger, and the ob- 
served direct continuum (i.e. powerlaw) stronger. The line profile is 
narrower than in regime I because the emissivity profile of the disk 
is flatter, as result of a nearly isotropic illumination. For NGC 5408 
X-1 and HOLM II X-I, we find a clear reflection component but 
the steepness of the spectra means that the Fe K features are less 
well-measured. 

The spectral solution we have found here for ULX suggests 
a new explanation for accretion onto spinning black holes. Our 
model assumes that the dominant source of radiation is a power- 
law continuum produced a few gravitational radii above a rapidly 
spinning black hole. Little thermal radiation is produced by the disc 
and is undetected by the current data. We assume that the power 
for this source is extracted magnetically from the disc and trans- 
ferred to the emission region by magnetic fields. Some of the power 
may even be extracted mag netically from the spin of the black hole 
dBIandford & Znaie kll977l) . Since radiation is only produced above 
the disc, radiation pressure need not oppose accretion. Indeed it will 
help squash the disc and maintain the high surface density required 
for our relatively low ionization parameters and thus observable re- 
flection. (Note that strong light bending occurs in the region being 
considered which is very close to a black hole, meaning that there 
is beaming but in the opposite sense to that normally envisaged.) 

The relevant radiation for computing the physical Eddington 
limit in this situation is the thermal disc radiation, not the power- 
law continuum and reflection associated with it. Since we detect no 
such thermal radiation, then the situation may be sub-Eddington, 
even for stellar mass black holes. Whether this solution can work 
depends on the extent to which magnetic energy extraction can be 
clean, in the sense of not requiring considerable thermal energy 
release. We note that the accretion flow is super-Eddington in the 
conventional interpretation in which the total energy release is con- 
sidered (especially since some radiation falls straight into the black 
hole). 

If the above solution is appropriate for these objects then they 
can either be stellar mass black holes with masses of kIOMq or 
IMBH of IOOsMq. Their spin is then likely to be natal. If they 
form through stellar collapse and the progenitor star was highly 
spinning, then it would be likely that the black hole would remain 
highly spinning. Many studies have argued that a binary compan- 
ion can spin up a massive sta r but the ma gnitude of the spin tip 
is still a matter of debate (e.g.. |Paczvnski e t al. 1998; Frver 1993; 
iBrown et alJl200Q ). Accretion would presumably come from a bi- 
nary companion. An interesting alternative is that they accrete from 
a debris disc pro duced by fallback from the (failed) supernova or 
coUapsar ( 101200 3). This could account for the high metallicity we 
infer (note that hydrogen-poor accretion would also slightly raise 
the Eddington limit). 

In summary, some ULX may consist of fast spinning black 
holes accreting rapidly from a metal-rich disc. Accretion power 
passes magnetically to the primary emission region which lies just 
outside the accretion flow. Strong reflection occurs, creating a soft 
excess and very broad iron-K features in the 0.3-10 keV X-ray 
spectrum. Sensitive hard X-ray spectra in the 10-50 keV band may 
test this scenario. 
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NGC 1313 X-1 NGC 1313 X-2 




Energy (keV) Energy (keV) 

Figure 2. The best-fitting (blurred) ionized reflection plus power-law model plotted on the spectra of NGC 1313 X-1, NGC 1313 X-2. M 81 X-6, HOLM IX 
X-1, NGC 5408 X-1 and HOLM II X-1 (upper-left to lower-right). Simple photoelectric absorption has also been applied. The model fits the data very well 
over the entire observed energy band. Data have been rebinned for visual clarity. 
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Table 2. Spectral parameters from fits with tiie laor and tlie reflection model. 



Parameter 


NGC 1313 X-1 


NGC 1313 X-2 


M 81 X-6 


HOLM IX X-1 


NGC 5408 X-1 


HOLM II X-1 


Laor Model 


E (keV) 


6.79±0.18 


6.87±0.10 


6.85±0.13 


6.91±0.06 


6.48±0.08 


6.4±0.7 


EW (keV) 


0.36 


0.45 


1.27 


0.70 


0.43 


0.15 


9 


7.1±0.6 


9.1±0.5 


9.95±0.05 


6.8±0.5 


9.77±0.23 


8.05±1.9 


T 


1.77 ±0.03 


2.02±0.05 


2.56±0.02 


1.57±0.02 


2.94±0.09 


2.61±0.03 




645/680 


709/711 


646/599 


1000/1036 


380/335 


578/539 


Blun'ed Reflection Model 




0.241±0.002 


0.374±0.003 


0.422±0.004 


0.164±0.002 


0.129±0.001 


0.179±0.003 


r 
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1.504±0.005 
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1.504±0.004 


3.04±0.01 


2.72±0.0I 
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4.61 ±0.07 


9.02±0.10 


8.72±0.09 


4.52±0.14 


6.5±0.6 


4.8±1.1 


Rin {Rg = GM/c2) 


1.36±0.07 


1.244±0.009 


1.314±0.008 


2.09±0.06 


2.45±0.08 


2.30±0.20 


Ro\it {Rg) 


400 


400 


400 


400 


400 


400 


iC) 


44.0±0.7 


46.2±0.4 


40.31±0.19 


50±10 


55.7±0.5 


59.5±0.5 


Fe/solar 


7.1±0.3 


7.7±0.5 


9.6±0.4 


4.98±0.12 


4.06±0.12 


4.5±0.8 


^{=47TF,,Jn) 


2155±12 


2531±13 


1789±14 


1673±19 


49.4±0.8 


29±6 


^powerlaw / ^total ^ 

i(0.05-ioo)keV (ergs"') 


0.32 
975/919 
1.4x10'*° 


_b 

991/952 
2.1xl0-*» 


956/840 
1.9x10"" 


0.77 
1292/1275 
3.6x10"" 


0.75 
736/577 
6.7x10"" 


0.92 
945/779 
1.4x10"" 


a(= cJIGM-) 


0.960 - 0.998 


0.970 - 0.998 


0.993 - 0.998 


0.95±0.17 


70+" '- ^ 

"■'"-0.05 — 


0.47±0.12£ 



" Fluxes calculated in the 0.05-100 keV energy band. 

^ For NGC 1313 X-2 and M 81 X-6 the best fits were obtained considering fully reflection dominated models, i.e. fpowcilaw/^total = 0- 
kerrconv gave a worse fit (755/578) that kdblur for this source, 
kerrconv provides a fit of the same quality (945/779) as kdblur (see text). 
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